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— . (57) Abstract: A method for determining tissue type includes quantitatively determining a tissue blood flow (TBD) deconvoluting 
2 Q(t) and C B (t), where Q(t) represents a curve of specific mass of contrast, and C 8 (t) represents an arterial curve of contrast concen- 
tration, and quantitatively determining a tissue blood volume (TBV) by deconvoluting Q(t) and C 8 (t). The method also includes 
quantitatively determining a tissue mean transit time (TMTT) by deconvoluting Q(t) and C,(t), and quantitatively determining a 
tissue capillary permeability surface area product (TPS) by deconvolving Q(t) and Q,(t). The method also includes determining a 
^ tissue type based on the TBF, the TBV the TMTT, and the TPS. 
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METHOD AND APPARATUS FOR CALCULATING 
BLOOD FLOW PARAMETERS 



CROSS REFERENCE TO RELATED APPLICATIONS 

» 

This application claims the benefit of U.S. provisional application No. 
60/243,196 filed October 25, 2000. 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to methods and apparatus for 
imaging systems and, more particularly, for calculating blood flow parameters. 

[0002] Stroke is a leading cause of disability among adults in North 
America. A National Institute for Neurological Diseases and Stroke (NINDS) trial 
was the first study to demonstrate that thrombolytic treatment within the first 3 hours 
of symptom onset is beneficial to stroke victims. However, treating patients beyond 
the initial 3 hours of illness with a thrombolytic treatment may carry a risk of 
intracranial bleeding which has resulted in a setback in the use of thrombolysis. This 
setback in the use of thrombolysis has prompted the realization that a more rigorous 
selection of patients according to certain criteria may reduce the risk of intracranial 
bleeding and thus, extend the therapeutic window for this therapy beyond the current 
3-hour limit to benefit more patients. 

BRIEF SUMMARY OF THE INVENTION 

[0003] In one aspect, a method for determining tissue type is 
provided. The method includes quantitatively determining a tissue blood flow (TBF) 
by deconvoluting Q(t) and C a (t), where Q(t) represents a curve of specific mass of 
contrast, and C a (t) represents an arterial curve of contrast concentration, and 
quantitatively determining a tissue blood volume (TBV) by deconvoluting Q(t) and 
C a (t). The method also includes quantitatively determining a tissue mean transit time 
(TMTT) by deconvoluting Q(t) and C a (t), and quantitatively determining a tissue 
capillary permeability surface area product (TPS) by deconvoluting Q(t) and C a ({)- 
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The method also includes determining a tissue type based on the TBF, the TBV, the 
TMTT, and the TPS. 

[0004] In another aspect, a system including at least one of a 
computed tomography system and a nuclear magnetic resonance system is provided. 
The system is configured to quantitatively determine a tissue blood flow (TBF) by 
deconvoluting Q(t) and C a (t), where Q(t) represents a curve of specific mass of 
contrast, and C a (t) represents an arterial curve of contrast concentration, and 
quantitatively determine a tissue blood volume (TBV) by deconvoluting Q(t) and 
C a (t). The system is also configured to quantitatively determine a tissue mean transit 
time (TMTT) by deconvoluting Q(t) and C a (t), and quantitatively determine a tissue 
capillary permeability surface area product (TPS) by deconvoluting Q(t) and C a (t). 
The system is also configured to determine a tissue type based on the TBF, the TBV, 
the TMTT, and the TPS. 

[0005] In another embodiment, a computer readable medium 
encoded with a program executable by a computer for processing scanning data is 
provided. The program is configured to instruct the computer to quantitatively 
determine a tissue blood flow (TBF) by deconvoluting Q(t) and C a (t), where Q(t) 
represents a curve of specific mass of contrast, and C a (t) represents an arterial curve of 
contrast concentration, and quantitatively determine a tissue blood volume (TBV) by 
deconvoluting Q(t) and C a (t). The program is also configured to instruct the computer 
to quantitatively determine a tissue mean transit time (TMTT) by deconvoluting Q(t) 
and C a (t), and quantitatively determine a tissue capillary permeability surface area 
product (TPS) by deconvoluting Q(t) and C a (t). The program is also configured to 
instruct the computer to determine a tissue type based on the TBF, the TBV, the 
TMTT, and the TPS. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0006] Figure 1 is a pictorial view of a CT imaging system. 
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[0007] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

[00OS] At present, there are both experimental and clinical data to 
suggest that the level of residual blood flow in the ischemic zone could be a useful 
indicator of the risk of intracranial bleeding. 

[0009] The importance of absolute measurement in stroke is also 
predicated by the fact that there exist distinct thresholds of cerebral blood flow (CBF) 
for various functions of the brain. It is possible to use these thresholds to decide 
whether a particular brain region is salvageable, hence thrombolytics should be used 
to restore CBF, or is already infarcted such that thrombolysis would achieve little but 
instead would increase the risk of intracranial bleeding. 

[0010] The brain has an intricate system of control to maintain 
cerebral normal limits when cerebral perfusion pressure decreases but remains within 
certain limits. This is achieved by dilation of the resistance blood vessels, i.e. the 
arterioles, vascular resistance and an increase of the cerebral blood volume (CBV), 
The concept of autoregulation of CBF can be used to advantage in stroke. For 
ischemic but viable tissue, autoregulation would lead to increased CBV so that mean 
transit time (MTT) is prolonged since MTT is the ratio of CBV and CBF because of 
the Central Volume Principle. On the other hand for ischemic but nonviable tissue, 
autoregulation is abolished such that both CBV and CBF are reduced but MTT may 
remain normal or only slightly elevated. 



[0011] In summary, the absolute measurements of CBF, CBV and 
MTT would permit the distinction between salvageable and infarcted tissue by the 
following criteria: 



Ischemic Tissue 
Type 


CBF 


CBV 


MTT 


Viable 




+ 


++ 


Infarcted 






-/+ 
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[0012] Wherein - designates a large decrease, - designates a slight 
decrease, + designates a slight increase, ++ designates a large increase, and -/+ 
designates either a slight increase or a slight decrease. 

[0013] It is the mismatch between CBF and CBV that discriminates 
salvageable and infarcted tissue. The corollary of this is that measurement of CBF 
alone will not reliably differentiate between viable and non-viable ischemic tissue. In 
addition, the ability to monitor quantitative changes in CBV may obviate the need to 
perform supplementary tests to evaluate the 'vascular' reserve of a stroke patient. In 
brief, the classical test of 'vascular' reserve evaluates the change in CBF in response 
to increases in pC02 (partial pressure of C02) in brain tissue which can be induced 
either by increasing inspired air concentration of C02 or by intravenous 
administration of a drug, such as a carbonic anhydrase inhibitor such as, for example, 
Diamox. For a normal brain, raised pC02 in tissue would induce a large increase in 
CBF. For ischemic tissue, where autoregulation is still intact, because CBV is already 
activated, the increase in CBF will be attenuated. Thus, tissue with a positive reserve 
is viable whereas that with little or no reserve left is at risk of infarction. 

[0014] As used herein Cerebral Blood Flow (CBF) is the volume 
flow of blood through vasculature including the large conductance vessels, arteries, 
arterioles, capillaries, venules, veins and sinuses. It typically has units of 
ml/min/lOOg. Also as used herein, Cerebral Blood Volume (CBV) is the volume of 
blood in the vasculature including the large conductance vessels, arteries, arterioles, 
capillaries, venules, veins and sinuses. It typically has units of ml/g. Additionally, as 
used herein Mean Transit Time (MTT) references that blood traverses the vasculature 
through different pathlengths such that there does not exist an unique transit time from 
the arterial inlet to the venous outlet. Instead there is a distribution of transit times 
and the mean transit time is the mean time of such a distribution. Furthermore, 
Minimum Transit Time (TTmin), as used herein, is the minimum time interval 
between entry at the arterial inlet and exit at the venous outlet of blood or contrast 
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medium. The Central Volume Principle relates the above three quantities in the 

CBV 

following relationship: CBF = . 

[0015] Referring to Figures 1 and 2, a multi-slice scanning imaging 
system, for example, computed tomography (CT) imaging system 10, is shown as 
including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a detector 
array 18 on the opposite side of gantry 12. Detector array 13 is formed by a plurality 
of detector rows (not shown) including a plurality of detector elements 20 which 
together sense the projected x-rays that pass through an object, such as a medical 
patient 22. Each detector element 20 produces an electrical signal that represents the 
intensity of an impinging x-ray beam and hence the attenuation of the beam as it 
passes through object or patient 22. During a scan to acquire x-ray projection data, 
gantry 12 and the components mounted thereon rotate about a center of rotation 24. 
Figure 2 shows only a single row of detector elements 20 (i.e., a detector row). 
However, a multislice detector array 1 8 includes a plurality of parallel detector rows 
of detector elements 20 so that projection data corresponding to a plurality of parallel 
slices are, or can be acquired simultaneously during a scan. 

[0016] Rotation of gantry 12 and the operation of x-ray source 14 are 
governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 
analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high-speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

[0017] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
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cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The qperator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, x-ray 
controller 28 and gantry motor controller 30. In addition, computer 36 operates a 
table motor controller 44 which controls a motorized table 46 to position patient 22 in 
gantry 12. Particularly, table 46 moves portions of patient 22 through gantry opening 
48. In one embodiment, computer 36 includes a device 50, for example, a floppy disk 
drive or CD-ROM drive, for reading instructions and/or data from a computer- 
readable medium 52, such as a floppy disk or CD-ROM. In another embodiment, 
computer 36 executes instructions stored in firmware (not shown). Computer 36 is 
programmed to perform functions described herein, accordingly, as used herein, the 
term computer is not limited to just those integrated circuits referred to in the art as 
computers, but broadly refers to computers, processors, microcontrollers, 
microcomputers, programmable logic controllers, application specific integrated 
circuits, and other programmable circuits. 

[0018] In one embodiment, system 10 is used to perform CT scans 
and determines blood flow parameters such as Tissue Blood Flow (TBF), Tissue 
Blood Volume (TBV), Tissue Mean Transit Time (TMTT) and Tissue Capillary 
Permeability Surface Area Product (TPS) as described below. In another 
embodiment, a Nuclear Magnetic Resonance (NMR) system (not shown) scans and 
determines TBF, TBV, TMTT and TPS blood flow parameters such as described 
below. In an exemplary embodiment, system 10 scans and determines cerebral blood 
flow parameters such as Cerebral Blood Flow (CBF), Cerebral Blood Volume (CBV), 
and Cerebral Mean Transit Time (CMTT). 

[0019] In one embodiment, system 10 is used to determine tissue 
type. More specifically, a tissue blood flow (TBF) is quantitatively determined by 
deconvoluting Q(t) and C a (t), where Q(t) represents the tissue residue function and is a 
curve of specific mass of contrast in tissue, and C a (t) represents an arterial curve of 
contrast concentration. Also a tissue blood volume (TBV), a tissue mean transit time 
(TMTT), and a tissue capillary permeability surface area product (TPS) are 



WO 02/079776 PCT/1B01/02898 



quantitatively determined by deconvoluting Q(t) and C a (t). In another embodiment, a 
Nuclear Magnetic Resonance (NMR) system (not shown) scans and determines tissue 
type by quantitatively determining TBF, TBV, TMTT and TPS. 

[0020] In one embodiment, the arterial curve of contrast 
concentration measured by system 10 and NMR system (not shown) is corrected for 
partial volume averaging as described herein. For example, during a cranial scan, but 
not limited to cranial scans, arterial regions within the vascular territories of the 
cerebral arteries (anterior and middle) are identified, and used to generate the 
measured arterial curve of contrast concentration, C a (t) . The measured arterial curve 
of contrast concentration is related to the arterial curve of contrast concentration, C a (t) 
byC fl , (0 = ^*C J (0> where k is the partial volume averaging scaling factor as 
explained in greater detail below. A venous region either within the sagittal or 
tranverse sinuses is located, and C v (t) is generated where C v (t) = C a (t)*h(t) and h(t) 

is the transit time spectrum of the brain, as explained herein. C a (t) and C v (t) are 
deconvolved to find ^p-. And a trailing slope of C a (t) is extrapolated with a 

m 

monoexponential function to find C aex (t) which is convolved with k to find 
C v ,«c(0- Where k is a partial volume averaging (PVA) scaling factor and is 

00 

\Ca,cx(t)dt 

determined according to k = \ . The measured arterial curve of contrast 

jC*,ev(0 dt 
0 

concentration, C' a (t) , is then corrected for partial volume averaging by dividing with 
the factor k to arrive at the arterial curve of contrast concentration, C a (t). 

[0021] A CBF and a CBV functional maps are generated with C? 0 (t) 
and divided by the PVA scaling factor, k. 

[0022] In one embodiment, a convolution integral for a tissue residue 
function Q(t) is: 
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Q(t)= JC a (t-u)h(u)du 

0 

[0023] where C a (t) is an arterial curve of contrast concentration curve 
and h(t) is an impulse residue function. Additionally, Q(t) is equal to the convolution 
of C a (t) and h(t) when the system is linear and stationary. Note that when C a (t). = 6(t - 
u), then the tissue residue function Q(t) is: 

t 

Q(t)= J8(t-u)h(u)du = h(t) 

0 

[0024] This relationship lends itself to the following interpretation of 
the convolution integral, e.g.C fl (f)= J^C a {t-nAt)5(t -nAt), therefore the tissue 

n 

residue function Q(t) is: 
t 

Q(t) = jC 0 (t- u)h(u)du , which is 

0 

t 

= l£[C.(t-u)8(t-u-nto)Xu)dit , which is 

0 « 

I 

= E f C « (' - ")<^ - " - nAt)h{u)du , which is 

" o 

[0025] A discretization of the tissue residue function convolution 
integral is determined by dividing a time interval [0,t] into m equal intervals of 
length At, each using the rectangular rule for quadrature of the convolution integral, 

according to: Q(t) = Q(mAt) = ^ C a ([ m - i]At)h(iAt)At . 



[0026] In a matrix notation, the convolution integral for the tissue 
residue function after discretization is Q = Ah and therefore Q = Ah is: 
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qM 




Q((M-l)^t) 




,q(mai) j 





c 8 (ai) 

C a (2A.) 



C a ((M-l)w) C a ((M-2)At) ... 

C a (MAl) C a ((M-l)Al) C a ((M-2)dl) 



C a (M-N)At) 
C a ((M-N + l))At 





fh(o) ^ 




h(At) 




h(2il) 




^h((N-l)Al) ; 



where Q is a M x 1 vector, A is a M x N matrix, and h is a N x 1 vector. A 
least squares solution, h , of Q = All facilitates minimizing the norm of a tissue 

residue vector, r = Q - Ah according to: ||r|| 2 = |q - Ahj 

[0027] In one embodiment, an equality constraint is incorporated into 
the least squares problem such that when h is the impulse residue function of a linear 
and stationary flow system, it will satisfy the time causality, i.e. if the injection of 
tracer does not occur at time zero, then some beginning elements of h should be set to 
zero (time causality), and defining a minimum transit time, i.e. h will have a plateau of 
a duration equal to the minimum transit time. The requirement for time causality and 
minimum transit time can be written as the following equalities: 

hj =o 



h, =0 

h I+2 = h l+3 



[002S] where the first I elements of h are zero, and the duration of the 
plateau is (L-l) dt, where dt is the sampling interval of each element of h. The 
equality constraints can be written compactly as the matrix equation Ch = b: 
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i 


N-l-L 




1 


0 •• 


0 


0 


0 


0 •■• 


0 0 








0 


1 •• 


0 


0 


0 


0 •■• 


0 0 








0 


0 •• 


1 


0 


0 


0 •• 


0 0 


o 


0 




0 


0 ! 


0 


1 


-1 


0 • 


•00 


o 


0 


1 


0 


0 : 


0 


0 


1 


-1 • 


•00 


o 


0 




0 


0 •• 


0 


0 


0 


0 • 


.. 1 _1 


o 


0 



where C is a Mix N matrix, h is the impulse residue function, and b=0 is a 
Nxl zero vector. Note that M)=I+L-1. 



[0029] Also C can be partitioned as: 



1 0 
0 1 

0 0 



0 0 0 
0 0 0 



0 
0 



10 0 0 



0 


0 


0 


0 


0 


0 


0 






0 


0 


0 


0 


0 


0 


0 


0 


0 



L-l 



0 0 
0 0 



0 1 

0 0 



-1 0 
1 -1 



0 0 



0 0 0 



1 -1 0 



where Ci is a Mi x Mi square matrix with full rank, and C2 is a M t x (N-Mi) matrix 
with only one nonzero element, wherein the nonzero element (Mi+1, l)th element is 
equal to -1 . Similarly h can be partitioned as: 



h = 



, where hi is a M) x 1 vector of the first Mi elements of h, and h2 is a N-M| x 
1 vector of the remaining elements. 

[0030] With these partitions of C and h, the constraint equation, Ch = 
b, can be written as: [c, j = b => c,ii, = (b - C,/; 3 ). Since Ci is square and of 
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full rank, its inverse exists, then the constraint equation can be written as: 
A, = C," 1 (b - C 2 h 2 ) , where the inverse of Ci is: 




L-l 



L-\ 



0 0 



1 



0 0 0 0 0 



0 0 0 0 
0 0 0 0 



0 0 0 0 
0 0 0 0 



0 0 0 0 0 
0 0 0 0 0 



1 1 

0 1 



1 1 
1 1 



0 0 0 1 1 
0 0 0 0 1 



[0031] As a check, hi = C^'(b - C2I12) is evaluated to ensure that it 
does not violate the constraint described herein. For example, since b is a zero vector, 
hi can be simplified as: hi =- C," 1 C2I12. As stated previously herein, C2 is a Mi x N-Mi 
matrix of the form: 



C,= 



0 


0 •• 


0 


0 


0 ■• 


0 


0 


0 •• 


0 


0 


0 •• 


0 


0 


0 • 


0 


0 


0 •• 


0 


-1 


0 ••• 


0 



1 




0 


0 


0 ••• 


0 




0 


0 


0 ••• 


0 




"0 0" 


0 


0 


0 ••• 


0 




-1 0 


-1 


0 


0 ••• 


0 





(h 



and I12 can be partitioned as: h 2 = 



\ ,l N J 
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A/.-l 



then,C 2 /i 2 = 



1 tf-A/, 

0 0 
-1 0 



V.1 


1 _ 


' 0 1 


, /J 2 > 


Af-A/,-1 





A/.-l 



and, 



lt l*2 



u 



=-c, 



/+£/ 



£-1 



A/,-1 0 

( 0 
1 0 



-1 
-1 



' 0 



A/,-1 



' o x 

0 
0 



as the time causality and minimum transit time constraints required. 

[0032] In one embodiment, the equality constraint Ch = b affect on 
the solution of the least squares problem Q = Ah can be determined according to: 



First, A is partitioned in the same way as C: A = A/ 

00 f 



A/, /V-A/, 

A x , Ay 



as: ,4/7 = [4,4] 



= AJ\ + A 2 h 2 = 4 



, then Ah can be written 



a 2 -.a } c-;c 2 \h 2 



Therefore with the equality constraint, Ch = b, the least squares problem All = Q is 
equivalent to the 'reduced' problem of (a, - AjC^C 2 )i 2 = Q or in short 

A r h 2 =Q where A r = (a 2 - A J C; , C 2 ). 

A 

[0033] From the least squares solution, 1)2, of the reduced problem 
A r h 2 = Q , the flill impulse residue function, h, can be reconstituted by : 
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h = 



v h 2/ 



and h, is a (I+L-l) x 1 vector whose first I elements are zero and the 



subsequent L-l elements are equal to the first element of I12 . 

[0034] For example, for the evaluation of A 2 - A } C~ l C 2 















1 


' X 


0 




0 








0 


M,-l 


0 


0 


= M,-1 


X 


0 
1 




0 








1 


1 


-1 


0 0 



= M,-1 





N-M,-) 




T 


0 0 




0 




0 


0 




-1 




L-l 


-1 







0 






0 


0 


0 




-1 






0 


-1 





where X is a plurality of partitions of matrices that are not used. Ai is then partitioned 
according to: 







L-l 


M-L + l 


X 


T . 


L-l 


X 


T 2 



such that: 
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A,C-'C 2 = 





M,-L+l 






1 


N-M,-l 














M-L+l 


X 


T, 


M|-L+l 


0 


0 










0 












-1 




L-l 


X 




L-l 


-1 


0 



M-L+l 



L-l 




therefore, ArA^C^C^ is the same as A 2 except that its first column is modified by 
adding the M x 1 vector 

rn 



M-L+l 



L-l 



to it. 

[0035] In one embodiment, a smoothness constraint by Lagrange 
Multiplier is incorporated into the least squares solution of the reduced problem 
A r h 2 = Q . In fact, if h is the impulse residue function of a physiologically realizable 
flow system, then its elements would change smoothly. In one embodiment, a matrix 
?f is appended to A r = (A 2 - AjC" 1 C 2 ) and vector yd to Q to the reduced problem to 
obtain: 



(a} 

h-> <=> p = Eh 0 where p = is a M2XI vector 

[yd J 

is a M 2 x (N-Mi) matrix, andp = Eh 2 is referred to the appended 





"A r " 







and E = 



A, 
yS 

reduced problem. 
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[0036] The least squares solution, h 2 , of the appended reduced 
problem p = Eh 2 facilitates minimizing the normal of the residual vector according 
to: 



2 

J J 



'Q 

yd 



yS 





"A," 






J 







=fc T 

= [Q T 

= Q T Q+T 2 d T d-hXQ- 7 2 hIS T d-Q T A r h 2 -y 2 d T Sh 2 '+fijAjA,fi a +y 2 hIS T Fh 
= Q T Q-hjA?Q-q T A r h, + fijA*A r h, +7 2 [d T d-h^F T d-d T Fh, + hjF T Fh,] 
Q-A r h 2 f +Y 2 ||d-Sh 2 



vS 



where y is the Lagrange multiplier and determines the relative weighting the least 



squares solution of the reduced problem places in minimizing |Q-A r h 2 | and 

more than 



||rf-F/? 2 | . If y is large, then the solution will minimize ^d-Fh 2 

[0037] The smoothness of the least squares estimate h 2 of the 
appended reduced problem can be gauged by the norm of its second derivative with 
respect to time according to: 







' l h (0) ^ 






4(2) 


dh 2 _ 1 
dt At 


h 2 (2)-h 2 (\) 


_ 1 

~ At 















4(1) 

h(2)-h(\) 
{h 2 (N-M i )-h 2 (N-M l 



-15- 



WO 02/079776 



PCT/1B01/02898 



d?~A? 



f^N-M^-lh^N-M^ -\)+ii£4-M x -2j 



"a? 



/i,^-M 1 ;-2^-Jwj-v+/i/N-M 1 -2; 



_j_ 

J/ 2 



10 0 0 
-2100 
1-210 
0 1-21 

0 0 0 0 
0 0 0 0 



0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 



0 1 -2 
0 0 1 



1 0 
-2 1 



} h (\) 

K(v 
kv 

f h (N-M x -2) 
h(N-M x ) 



To incorporate the smoothness constraint, the matrix S is set according to: 



S = - 



At 2 



10 0 0 
-2100 
1-210 
0 1-21 



0 
0 



0 0 0 
0 0 0 



0 0 
0 0 
0 0 



0 
0 
0 



0 0 0 

0 1 -2 
0 0 1 



0 0 

0 0 

0 0 

0 0 

1 0 
-2 1 



where d=0. The least square solution of the appended reduced problem facilitates 



minimizing both Q - A r h J and 



Sh, 



with the relative weighting controlled by y . 



[0038] In one embodiment, an inequality constraint is imposed on the 
solution of the appended reduced problem. The appended reduced problem is: 



r q N 




X" 


J d > 







which can be written as: p = Eh 2 . 

[0039] In another embodiment, the least squares estimate of the 
appended problem under the linear inequality constraint Gh 2 is determined. The 
dimension of the appended matrix E is M 2 x(N-M,). In one embodiment, E has a 
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rank of K<(N-M,) and has the following singular value decomposition: 



E = U 



S 0 
0 0 



V T where S is a K x K diagonal matrix with its diagonal entries equal to 



the singular values of E. U and V are orthogonal matrices of dimensions M 2 xM 2 and 
(N-M,)x(N-M,) respectively. The least squares solution, h 2J of the appended 
reduced problem p = Eh 2 will minimize ||Eh, -p||" . 



[0040] Using the singular value decomposition of E, |Eh 2 -p| 2 can 



be simplified as 
|Eh 2 -p|| 2 =U 



t 



S 0 
0 0 

S 0 
0 0 



H2 J 


U T U 

0 0 


V T h 2 -U T p 




"S 0" 


"v, T " 






2 




0 0_ 


_ V 2_ 




P 



SV, T h 2 
0 



u; P 
uTpJ 



SV, T h 2 -U^p' 

-uTp 



HKh.-u^ir+iuipi 



[0041] Since |ujp| 2 is constant, minimizing ||Eh 2 -p| 2 is 
equivalent to minimizing |SV 1 T h 2 - U^pll . 



[0042] In one embodiment, the variable Yi 2 is changed to z such 
thatz = SV 1 T h 2 -Ufp <=> V, T h 2 =S- 1 z + S- | U[p . 

[0043] Minimizing |Eh 3 -pf is the same as minimizing ||z|| 2 and, 
therefore the inequality constraint becomes: 
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Gh 2 >b 
O G(V,V, T +V 2 V 2 T )h 2 >b 
<=> GV.CS-'z + S-'U^ + GVjVjh, >b 
o GV,(S-'z + S-'U^p) + GV, V 2 T V,(S"'z + S _, U*p) > b 
o GV,S-'z+ GV 2 V, T S-'z > b - GV,S-'U^p - G^V/V^'U^p 
<=> G(I 2 +V 2 V 2 T )V,S- , z>b-G(I 2 +V,V 2 T )V 1 S- 1 U^p 

[0044] Therefore, the problem of minimizing |Eh 2 =pf subject to 
the constraint Gh 2 >b is equivalent to minimizing ||r|| 2 subject to 

Gfli+VjV^V.S-'z^b-Gflj+VjV^V.S-'Ufp. From the solution, z, that 

minimize |jz|| 2 , the corresponding least squares solution, h: , of the appended reduced 
problem Eh 2 = p can be found by the inversion of the equation z = SV, T h 2 -U|p . 

A 

From the least squares solution, h 2 , of the appended reduced problem, Eh 3 = p , the 
foil impulse residue function, h, can be reconstituted by : 



h = 



V ll2y 



and h, is a (I+L-l) x 1 vector whose first I elements are zero and the 



subsequent L-l elements are equal to the first element of I12 . 

[0045] In another embodiment, E has a rank of K=(N-Mi) and has 



the following singular value decomposition: E = U 



V T , where S is a K x K 



diagonal matrix with its diagonal entries equal to the singular values of E. U and V are 
orthogonal matrices of dimensions M 2 xM 2 and (N-M, ) x(N-Mj) respectively. The 

least squares solution, h 2 , of the appended reduced problem f = Eh 2 facilitates 

minimizing |Eh 2 -p|| 2 . 



[0046] Using a singular value decomposition of E, |Eh 2 -p|| 2 can 
be simplified as: 
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Eh 2 -p" - U 



V T h 2 -p| 



U T U 



V T h 2 -U T p| .Since U T is orthogonal 

2 





"s" 


II- 










2 ii 


"SV T h 2 " 










V T h 2 -U T p = 




V T h 2 - 




p 












0^ 










| 


o 




_UjpJ 



SV T h 2 -U^p 

-uJp . 



= SV T h 2 -U[p + ui 



since IjuTi:!' 2 is constant, minimizing|Eh 2 — p|j 2 is equivalent to minimizing 

||sV T h 2 — U7p| 2 • 1)1 one embodiment., the variable h? is chnnecd to z such that: 
z = SV T h 2 -Li; r p <=> V T h 2 = S- ] z 4- S' 1 Ujp 

[0047] Minimizing |jEh 2 -p|| 2 is the same as minimizing |z|| 2 and the 
inequality constraint becomes: 

Gh 2 >b 
o G(VV T )h 2 >b 

<=> GVS~'z>b-GVS- l Ufp 

[0043] Minimizing |Eh 2 -p|| 2 subject to the constraint Gh 2 >b is 
equivalent to minimizing ||z|| 2 subject to GVS _1 z > b -GVS"'U7"p . 

[00-!°] A formulation of smoothness, nonncp r;ive and monotonicity 
constraints fr- the reduced problem can now be stated. The h-ast squares solution of 
the reduced problem: Q = A r h 2 with the smoothness, nonnegntive and monotonicity 
constraints can be recast as the least squares solution of the appended reduced 
problem according to: 



q 

[yd 



ih, nr p = Eh, 



with constraint 3 (N . Mi)K(N . Mi) h 2 > 0 and, ])h 2 >0 } where D is a 



(N-M, -l)x (N-M,) matrix of ;he form: 
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D = 



1 -1 0 
0 1 -1 



0 0 0 
0 0 0 



0 0 0 
0 0 0 



1 -1 0 
0 1 -1 



such thatDh, s 0 is equivalent to h 2 (\) > h 2 (2 )>■■■ h 2 (N - MJ . 

The two inequality constraints I(n. M{)x(N . M Jh ^ 0 and Dh 2 so can be combined into 



h 2 >0 or, G/j 2 >0where G is a 2(N-M, -l)x(N-M,) matrix as 

'<N-M,-l).x(N-M0 " 

defined herein such that I (N _ Ml)x(N _ Mi) is an (N-Mi)x (N-Mi) identity matrix. 



[0050] The least squares solution of the appended reduced problem 
p = Eh 2 subject to the nonnegative and montonicity constraints can be determined 

using the method described herein. From the least squares solution, h2, of the 
appended reduced problem, Eh, =p, the full impulse residue function, h, can be 
reconstituted by : 



h = 



and h, is a (I+L-l) x 1 vector whose first I elements are zero and the 



subsequent L-l elements are equal to the first element of I12 . 

[0051] In one embodiment, in tissue where there is no leakage of 
contrast from the blood stream into the interstitial space, for example, the brain with 
intact blood-brain-barrier, the reconstituted full impulse residue function h(t) can be 
used to determine Tissue Blood Flow (TBF) as the peak height of h(t), Tissue Blood 
Volume (TBV) as the area underneath h(t) and Tissue Mean Transit Time (TMTT) as 
the TBV divided by TBF. 
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[0052] In one ' embodiment, in tissue where there is leakage of 
contrast from the blood stream into the interstitial space, the impulse residue function 
h(t) can be modeled by the Johnson- Wilson model. 

[0053] An adiabatic approximation of the Johnson- Wilson Model is 
then used to parameterize the impulse residue function in terms of TBF, TBV, TMTT 
and TPS. In one embodiment, the impulse residue function for the adiabatic 
approximation of the Johnson- Wilson model is: 



R(t) = 



0.0 t<t 0 

1.0 t 0 <t<t 0 +W (1) 

Ee -kO-o-w> t >t +W 



where t 0 is a time delay between an arterial curve of contrast concentration, C a (t), 
and the tissue residue function, Q(t), W is the tissue mean transit time (TMTT), and k 
is the rate constant of the exponential function and is defined according to: 

k=f (2) 

e 

where F is the Tissue Blood Flow (TBF), E is the extraction efficiency of 
contrast material from blood and V e is the distribution space of contrast material in the 
intersitial space. Tissue Blood Volume is equal to TBF x TMTT. Tissue capillary 
permeability surface area product (TPS) is: - ln(l - E) . 

[0054] The tissue residue function, Q(t), is given by: 

Q(t) = F*[C a (t)*R(t)] (3) 

where * is the convolution operation. 

[0055] The linearization of the tissue residue function (or the 
expression of the tissue residue function or some t function of it in terms of linear 
functions of the parameters (or some combination of them) of the model) is based on 
finding the (time) integral of Q(t) according to: 
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T t 



\Q{t)dt=\dt\C a {uW-u)du (4a) 

0 0 0 

[0056] Interchanging the order of integration: 



T-u 



jQ(t)dt = Jdu f C a (u)R(t - u)dt = Jc, (u)du |R(t - u)dt = Jc a (u)du {R(t)dt 

0 0 u 0 u 0 0 

(4b) 

T 

[0057] The time integral jQ(t>dt is evaluated for the three time 

0 

intervals 

,...x 0<T<t 
(in) (i) u - A - l o 

r , t Q <T <t D +W 

(iv) (n) 0 0 

(v) t 0 + W<T 

T 

[0058] For example when the time integral Jcxodt is evaluated for the 

for 

/in 

0<a<To0>-«>-T«T>T-ii>0«0<T-«<T<«., therefore: 



time intervalO<T<t 0 , then Je ( ,)</, = jc>)rf/jf*<0<* 



[0059] \Q{t)dt = Q (5) 

0 

T 

[0060] When the time integral jQ(t)dt is evaluated for the time 

0 

interval t 0 < T < t 0 + W , then setting T = t 0 +a 9 where a<W , the time integral is: 

T t 0 +a te+a-u 

\Q{t)dt = F \C„(u)du \R[t)dt 
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JC„ (ti)du + |C„ (w)rf« J/?(0 A = F. \C 0 (u)du J R(t)dt + F. ' JC <7 (u)rfu " Jrt(i)A 

o n Jv ° y j j o 



For 



0Zu<a,t B + a>t B +a-u>t o 0X t 0 + W >t o + a-u>t 0 , 



'then: 



Fjc»rf« \R{t)dt = F\CSti)du 

0 0 0 

= F\C a (u)du \R(t)dt = F(a-u)\C t (u)du 

0 /„ 0 

F* J C fl («)tfi< Jfl(f)<ft = 0 



[006 1 ] Combining A and B above: 



\Q{t)dt = F* \{a-u)C n (u)du or, 



(6a) 



(T-0 Jc»rf«- J«C»rf« =Jfi(0* 

o o Jo 

T 

[0062] When the time integral jQ(t)dt is evaluated for the time 



interval T > t 0 + W , then and setting T = t 0 + W + a , where a >0 then 



T-u 



fQ(t)dt = F.Jc,(u)du |R(t)dt 



= F* \C a (u)du \R{t)dt 

0 0 



= F 



t 0 +a+W 



Jc»du + Jc,(u)du+ j C,(u)du jR(t)dt 



a+W 
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=F| ) G(n)du ]R{t)dt + F [ C*{ii)du ]R{t)dt + F[^ C a (u)du [R(t)dt 



C: For O<u^a >to+ a + W>t 0 ^a + W-u>t c+ W then 



t„+a+IV-u 



FjC a (u)du I R(t)dt = FjC a (u)du 



t.+a+W-u 



jR(t)dt + J R{t)dt + J^(0^ 



= FPF Jc„ (it)du + F JC n («)rf« J* (0* 



= \C„ (u)du + F jC(u)du j Ee~ k{ '~'°~ W) dt = FW jc„ (u)du + FE Jc„ (n)du je'^dt 



= Fw)c a (u)du - )c a (u)du[e-^ - 1] 



= FW jc 8 («)rfw + — |C„ (u)du -— |c n (i/)e- t( "- u) rf» 

0 ^ 0 ^ C 0 



F: 



0(T) = F fC a (u)R(T - zO^/ = F \C a (u)R(t 0 + * + ff' - 



/,+tf+W 



= F JC fl («)jR(/, + a + FT - it)du + F \C a (u)R(t c + a + FF - + F Jc a + a + FF 



= FF JC fl {u)e' kia ' a) dit + F JC fl (u)rfu 

0 a 

Therefore, — fc„ (u)e^ (a ' u) du = - - \c a (u)du 
k x k k r 



and 



t_+a+W-u 



FjC a (u)du |R(t)dt = FW}C a (u)du + ^}C a (u)du-^+^ jC a (u)du 



Q(T) F : 



C: 0 



k k 



-24- 



WO 02/079776 



PCT/1B01/02898 



D: For 



a£u<a + W,t 0 + W-iiZt e +a + W-u>t 0 e>t 0 +WZt 0 +a + W-iiZt o 



F jC a (u)du j R{t)dt = F jC 0 (u)du 



(,+a+W-u 



fR(t)dt + \R{t)dt 



= F Jc„ (u)du J Jl(0 A = F[a + W - u] J C„ («i)dw 



E . For a + W £ u < a + W + t 01 t 0 Z t 0 + a + W - u > 0 

1,+a+IV l 0 +a+IV-u 

F jC a (tt)du jR(t)d( = 0 

a+W 0 

Therefore, for T - 1 0 + W + a , where a > 0 : 



= Ff'F \C a (u)du + — JC fl (w)tfu - MLi + _ J Cfl («)rfa + F[a + W-u] \C a (u)du 



(7a) 
or, 

01 = F/c 
= Fk 



0 T-r.-IP T-/ e -IK 



o 0 



T-t 0 -lV 



T-r e -/F 



i 

C a (u)du-kjQ(t)dt 



(7b) 



[0063] In one embodiment, the time integral { fi <<)rf ' f or the three 



time intervals 



(i) 



0<T</ 



(ii) t.STSt. + W 
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(iii) 



/ + W <> T 



generates the following system of linear equations 



-k\0{t)dt 
0= o 



0 = Fk 



0<T<t„ 



T-t 0 



(T-t 0 ) Jc.(u)du- JuC a (u)du 



Q(T) » Fk 



T-t 0 



T-< 0 



W | C a (u)du + (T-t 0 -w) | C a (u)du- juC : 



(u)du 



T-i 0 -W T-i 0 -W 



T-i 0 T-i 0 -W T 

-F Jc o (u)du + FE Jc a (u)du-k|Q(l)dt 



involving a plurality of unknowns Fk, F, FE and k, provided that to and W are known. 

[0064] For a set of given t 0 and W, a least squares algorithm is used 
to estimate Fk, F, FE and k subject to the following linear constraints: 

Fk>0 
F>0 
FE>0 
k>l 
Fk > F 
F > FE 

[0065] The algorithm for estimating to, W, Fk, F, FE and k includes 
two main golden section search subroutines. An outer golden section search 
subroutine searches for t 0 , an inner golden section search subroutine searches for W 
with the value of t 0 , assumed in the outer routine. Within the inner golden section 
search routine, values of t 0 and W are fixed and the optimization of Fk, F, FE and k 
can proceed with LDP algorithms. 

[0066] In one embodiment, values of F, E, Ve and k are: F=0.1 
ml/min/g, E=0.5, Ve0.25 ? interstitial space, and k=FE/Ve=0.2 min-l=12s-l. In 

-26- 



WO 02/079776 PCT/IBO 1/02898 



another embodiment, values of F, E, V e and k are: F=0.1 ml/min/g, E=0.25, V e 0.25, 
interstitial space, and k=FE/V e =0.1 min" l =6s' 1 . 

[0067] Although k>=l is the constraint, after solution, for any k > 5s" 
1 assumes that there is no leakage and do not calculate PS, and recalculate blood 
volume and MTT instead of taking W as MTT and EW as volume. The volume can 
be taken as the area underneath the flow scaled impulse residue function (F.R(t)), and 
MTT as the volume divided by the height of the flow scaled impulse residue function. 

[006S] Discretization of the system of linear equations is performed 

according to: 

[0069] Let At be the time interval, and t 0 =mAt, W=nAt, and T=iAt, 

then: 



/ 

k^QUW 0<i<m 

—in i-m j 

i-m-fi i-flm i 

+F ^Q0)^+FE^<^0)^-^Q[i)^ i^n+m 



o = n 



(/-m)A^C fl (;)A/-£ jC n (j)M 
/-I y=i 



QCI)=F1 



i-ni i-ni i-m 

j=4 j=i-in-n j=i-m-n 

(9) 



[0070] The various sums involving the arterial function Ca(t) can be 
pre-calculated and stored in 2-dimensional arrays AAR and AFM. For example, 

i-m 

^c 0 o>tas the [(i-m-n),(i-m)]th element of Area array referred to herein as an 

j=i-m-n 

i-m 

arterial area array (AAR). Also, ^]c, 7 (y)Aras the [(i-m-n),(i-m)]th element of First- 

jmi-m-u 

moment array referred to herein as an arterial first moment array (AFM). The various 
sums involving the tissue residue function Q(t) can be pre-calculated and stored in a 
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1 -dimensional array AQ. For example, ^0(/)a; as the ith element of an array for area 

jml 

of tissue residue function Q(t), (AQ). 

[0071] A solution for F, E, and k when to and W are known can be 
determined as follows: 



Using the arrays: Q, AAR, AFM and AQ, Eq. (9) can be written as Eq. (10): 



'0(0 = 0 




"0 




0 


0 


-AQ(\) 


«#/!)- 0 




0 




0 


0 


-AQ{m) 


Q(m + \) = 0 




bt.AAR(\>\)-AFM(\ t \) 




0 


0 


-AQ(m + \) 


g(m + 2) = 0 




TAt.AAR&l)- AFM(\ t 2) 




0 


0 


-AQ(m + 2) 


Q(M + n) 




nAt.AAR(ln)- AFM(\,n) 




0 


0 


-AQ(m + n) 


Q(m + n + )) 




nA.AAR(\,n + 1) + &i.AAR(\,n + 1) - AFM(\,n + 1) 




AAR(\ t \) 


-AQ(m + n + \) 


Q(m + n + 2) 




nAt.AARiln + 2) + 2AT.AAR{2,n + 2) - 


AFM(2,ii + 2) /M*(2,/i + 2) /M*(l,2) 


-AQ(m + n + 2) 


0(m + /i + p) 






- AFM(p,n + 


p)AAR(p i n 


+ p)AAR(\,p) 


-AQ(m + n + p) 



or Q = Ax 
and: 



QG)=0 




*0 


0 


0 


-AQ© 




Q(m)=0 




0 


0 


0 


-AQ(m) 




Q(m+1)=0 




At. A AR(,1) - AFM^l) 


0 


0 


-AQ(mfl) 




Q(m+2)=0 




2At.AAR(,2)-AFM(,2) 


0 


0 


-AQ(m+2) 




Q(M+n) 




nA.At.A>y,n)-AFjVJ{,n) 


0 


0 


-AQ(rofn) 


x = 


Q(m+n + l) 




n A. A. AAl,n+ 1) +At.AAR(,n + 1)- AFM$,n + 1) AAR(,n + 1) AARU) 


-AQ(m+n + l) 




Q(m+n + 2) 




n A.At. A^,n + 2) + 2.A AT. AAfyi + 2) - 


- A FMg.n + 2) A A R&n + 2) AAR{,2) - AQ(m+ n + 2) 




_Q(m+n+p) 




_nA.At.AAn + p)+ p A. At AARi + p) - AFM(p,m- p) AAR(p,rtf- p) AAR(,p)- AQ(m+ n + p) 





The constraints are: 
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1 


0 


0 


0 






0 


0 


1 


0 


0 


~Fk^ 




0 


0 


0 


1 


0 


F 




0 


0 


0 


0 


1 


FE 




1 


1 


-1 


0 


0 


_k 




0 


0 


1 


-1 


0 






0 



(11a) 



or, Gx £ b 



G = 



1 0 


0 


0 


0 1 


0 


0 


0 0 


1 


0 


0 0 


0 


1 


1 -1 


0 


0 


0 1 


-1 


0 



and 



b = 



[0072] Therefore, the solution for the model parameters given t 0 and 
W can be stated as linear least squares problem for x according to q = Ax, subject to 
the linear inequality constraints Gx >b wherein q is a (m+n+p) x 1 vector, A is a 
(m+n+p) x 4 matrix, b is a 6x1 vector, and x is a 4 x 1 vector of the model parameters. 
In one embodiment, A has a rank of 4, i.e., full column rank and has the following 



singular value decomposition: A = U 



0 



V T , where S is a 4 x 4 diagonal matrix with 



its diagonal entries equal to the singular values of A. U and V are orthogonal matrices 
of dimensions (m+n+p) x (m+n+p) and 4x4 respectively. 

[0073] The least squares solution, x , of the least squares problem Q 
= Ax will minimize ||Ax - Q||" . Using the singular value decomposition of A, 

|Ax-Q|"can be simplified as 



|Ax-Q|r=|U 



V T x-Q =U T U 



V T x-U T Q , since U is orthogonal, 
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V T x-U T Q = 



SV T x-U^Q 
-U*Q 



SV T x 
0 



L U 2Qj 



, and 



= ||SV T x-U^Q|r+||U^Q| 



[0074] Since ||uT Q|P ^ s constant, minimizing ||Ax-Q|| 2 is equivalent 

to minimizing ||SV T x -U^q|"*. In one embodiment, variable x is changed to z, 
therefore, 

z = SV T x-U 1 iQo V T x =S" , z + S" , U^Q. Minimizing |Ax-Qfis the same as 
minimizing ||z|| 2 . Therefore, the inequality constraint becomes: 

Gx>b 
«G((VV T )x>b 
oGv(S" l z + S" l U 1 [Q)*b 
oGv(s" , z + S-'U^Q)^b 
oGVS H z>b-GVS M U^Q 

[0075] Therefore, minimizing ||Ax-Q|| 2 subject to the constraint 

GX > b is equivalent to minimizing |z|" subject to GVS" ! z > b - GVS^U^Q . In one 

embodiment, the coefficient matrix of the above system of linear systems will be 
singular value decomposed (SVD). It can be seen that the coefficient matrix is 
dependent on Q(t), this means that the SVD of the coefficient matrix has to be 
repeated for each tissue residue function. The coefficient matrix is N x 4 where N is 
the number of data points of the tissue residue function. The first 3 columns of each 
coefficient matrix is the same being only dependent on the arterial function. 
Therefore, SVD of the N x 3 matrix can be calculated , and for each tissue residue 
function, the fourth column is appended according to: 



-*£2(y)A/,/c£e(y)^ 

;=i j=\ ;=i y«j 
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[0076] The SVD of the N x 4 coeffficient matrix is then calculated 
from an update of the SVD of the basic N x 3 matrix. In one embodiment, the effect 
of partial volume averaging (PVA) of the arterial curve on the determination of F and 
V with the deconvolution method can be determined. For example, letting Q(t) be the 
brain curve of specific mass of contrast in tissue, C a (t) be the arterial curve of contrast 
concentration, R(t) be the impulse residue function, F be the blood flow, V be the 
blood volume, and MTT be the mean transit time. Then, if blood flow is stationary 
and CT measurement is linear with respect to contrast concentration, by principle of 
linear superposition then Q(t) =FC a (t)*R(t). 

[0077] Also, if the arterial curve is underestimated due to partial 
volume averaging from the finite resolution of CT scanners (~ 8 lp/cm) then 
C„(t) = k.C a (t), where C a (t) is the arterial curve measured with partial volume 
averaging (PVA) and k is the multiplicative factor due to PVA, which is less than one. 

F . 

Therefore, Q(t) = —CJt)*R(t). The deconvolution between Q(t) and CJt) gives 
k 

F F F F V 

—R(t) . The maximum height of —R(t) is — and the area underneath — R(t) is — . 
k k * * k k 

This explains why when a partial volume averaged arterial curve is used in the 

deconvolution with the brain curve, both F and V are scaled by — ( >\) . Note that 

k 

F 

MTT = — is not scaled even in this case when the arterial curve is measured with 
partial volume averaging. 

[0078] In one embodiment, a correction for partial volume averaging 
of the arterial curve uses the factor k such that Cjt) = k* CJt), therefore 

OO 

\c'jt)dt 

] c = JL m ]sj 0 te that CJt) is the measured (known) arterial curve i.e. the integral 

\C a (t)dt 

0 

in the numerator can be calculated, while C a (t) is the true arterial curve which is 
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unknown. If we assume that there exist regions in veins which are free of partial 
volume averaging, then unlike arterial curves, venous curves can be measured 
accurately. 

[0079] For example, let h(t) be the transit time spectrum from artery 
to veins, therefore C v (t) = C a (t)*h(t) . Since h(t) is a transit time spectrum, by 

00 

definition jh(t)dt = l. Therefore jc/t)dt= \C a (t)dt , and * = , where 

0 0 0 \C v (t)dt 

o 

k is expressed in terms of the partial volume averaged arterial curve (Cj(t)) and the 
venous curve C v (t). Both of which are known (measured). In one embodiment, C a (i) 
and C v (t) do not return to baseline during the measurement time, therefore integration 

to infinity is not desirable since C v (t) = C (t) * h(t), and therefore, C v (t)=Cj't)*^. 

k 

In other words, ^^can be obtained by deconvolution between C v (t) zndC^O, for 
example if C (t) is the extrapolated arterial curve obtained by extrapolation of the 

a tar 

trailing slope of Cj'l) with an exponential function. C (t) returns to baseline very 
rapidly such that ^C at Jt)dt can be easily calculated. Since — — and C^Jt)m both 

known, their convolution can be calculated as C ve /# = C' (t)*^&. Similar to 

k 

C a ex (t),C v Jt) also returns rapidly to zero baseline and ^C v Jt)dt = — JC fl Jt)dt , or 

oo T 

\c'„Jt)dt \c' a Jt)dt 

/c = -2j = -^r , where [0,T] is the time interval within which 

\C y Jt)dt \C v Jt)dt 

0 0 

C^JOmd C V) ex(0 first increase from zero baseline and then decrease back to zero 
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baseline again. Therefore, the time shift between C' a Jt) and C v , cx (t) will not affect 
the calculation of k. 

[0080] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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WHAT IS CLAIMED IS: 

1 . A method for determining tissue type, said method comprising: 

quantitatively determining a tissue blood flow (TBF) by deconvoluting 
Q(t) and C a (t), where Q(t) represents a curve of specific mass of contrast, and C a (t) 
represents an arterial curve of contrast concentration; 

quantitatively determining a tissue blood volume (TBV) by 
deconvoluting Q(t) and C a (t); 

quantitatively determining a tissue mean transit time (TMTT) by 
deconvoluting Q(t) and C a (t); 

quantitatively determining a tissue capillary permeability surface area 
product (TPS) by deconvoluting Q(t) and C a (t); and 

determining a tissue type based on the TBF, the TBV, the TMTT, and 

the TPS. 

2. A method in accordance with Claim 1 wherein quantitatively 
determining a TBV comprises quantitatively determining a TBV for a tissue having a 
blood stream containing a contrast without leaking the contrast into an interstitial 
space of the tissue by solving a matrix equation of Q = All for a vector h, wherein 
vector h includes a plurality of elements comprising an impulse residue function at 
different times, Q comprises a vector including elements comprising values of a tissue 
residue function at different times, A comprises a matrix formed by values of the 
arterial curve of contrast concentration at different times, quantitatively determining a 
TBF comprises quantitatively determining a TBF for the tissue by solving the matrix 
equation of Q = Ah for the vector h, quantitatively determining a TMTT comprises 
quantitatively determining a TMTT for the tissue by solving the matrix equation of. 
Q = Ah for the vector h. 
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3. A method in accordance with Claim 2 further comprising 
determining a least squares solution for the vector h under an equality constraint. 

4. A method in accordance with Claim 3 wherein determining a 
least squares solution comprises determining a least squares solution for the vector h 
under a time causality constraint and a minimum transit time constraint. 

5. A method in accordance with Claim 4 wherein determining a 
least squares solution further comprises determining a least squares solution for the 
vector h under a smoothness constraint, a nonnegativity constraint, and a monotonicity 
constraint, wherein the smoothness constraint forces h to be smoothly varying, and the 
monotonicity and nonnegativity constraints force h to start at a maximum and then 
monotonically decreases towards a zero baseline. 

6. A method in accordance with Claim 1 wherein quantitatively 
determining a TBV comprises quantitatively determining a TBV for a tissue having a 
blood stream containing a contrast with leaking the contrast into an interstitial space 
of the tissue by solving a matrix equation of Q = Ax from the linearization of the 
tissue residue function for a vector x, wherein vector x includes a plurality of elements 
comprising TBF, TBV, TMTT ? TPS and combinations thereof, Q comprises a vector 
including elements comprising values of a tissue residue function at different times, A 
comprises a matrix formed by values of the arterial curve of contrast concentration 
and tissue residue function at different times and combinations thereof, quantitatively 
determining a TBF comprises quantitatively determining a TBF for the tissue by 
solving the matrix equation of Q = Ax for the vector x, quantitatively determining a 
TMTT comprises quantitatively determining a TMTT for the tissue by solving the 
matrix equation of Q = Axfor the vector x, quantitatively determining a TPS 
comprises quantitatively determining a TPS for the tissue by solving the matrix 
equation of Q = Ax for the vector x. 

7. A method in accordance with Claim 6 further comprising 
determining a least squares solution for the TBF, the TBV, the TMTT, and the TPS 
under a nonnegativity constraint that determines a TBV, a TBF, a TMTT and a TPS 
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for the tissue where there is leakage of the contrast from the blood stream into the 
interstitial space. 

8. A method in accordance with Claim 1 further comprising 
quantitatively determining a partial volume averaging scaling factor for the arterial 
curve of contrast concentration by: 

deconvolving the measured arterial curve of contrast concentration 
with a venous curve of contrast concentration to determine a transit time spectrum 
through a tissue of interest; 

extrapolating the arterial curve; and 

convolving the extrapolated arterial curve with the transit time 
spectrum to generate an extrapolated venous curve, wherein the partial volume 
averaging scaling factor is the ratio of an area underneath the extrapolated arterial 
curve to an area underneath the extrapolated venous curve. 

9. A method in accordance with Claim 1 further comprising 
scanning with at least one of a computed tomography (CT) system and Nuclear 
Magnetic Resonance system (NMR) to measure the arterial curve of contrast 
concentration and a tissue residue function, said quantitatively determining a TBV 
comprises quantitatively detemiining a TBV for a tissue having a blood stream 
containing a contrast without leaking the contrast into an interstitial space of the tissue 
by solving a matrix equation of Q = All for a vector h, wherein vector h includes a 
plurality of elements comprising an impulse residue function at different times, Q 
comprises a vector including elements comprising values of the tissue residue 
function at different times, A comprises a matrix formed by values of the arterial 
curve of contrast concentration at different times, said quantitatively determining a 
TBF comprises quantitatively determining a TBF for the tissue by solving the matrix 
equation of Q = All for the vector h, said quantitatively determining a TMTT 
comprises quantitatively determining a TMTT for the tissue by solving the matrix 
equation of Q = All for the vector h. 
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10. A method in accordance with Claim 1 further comprising 
scanning with at least one of a computed tomography (CT) system and Nuclear 
Magnetic Resonance system (NMR) to measure the arterial curve of contrast 
concentration and a tissue residue function, said quantitatively determining a TBV 
comprises quantitatively determining a TBV for a tissue having a blood stream 
containing a contrast with leaking the contrast into an interstitial space of the tissue by 
solving a matrix equation of Q = Ax from the linearization of the tissue residue 
function for a vector x ; wherein vector x includes a plurality of elements comprising 
TBF, TBV, TMTT, TPS and combinations thereof, Q comprises a vector including 
elements comprising values of a tissue residue function at different times, A 
comprises a matrix formed by values of the arterial curve of contrast concentration 
and tissue residue funnction at different times and combinations thereof, said 
quantitatively determining a TBF comprises quantitatively determining a TBF for the 
tissue by solving the matrix equation of Q = Ax for the vector x, said quantitatively 
determining a TMTT comprises quantitatively determining a TMTT for the tissue by 
solving the matrix equation of Q = Axfor the vector x, said quantitatively 
determining a TPS comprises quantitatively determining a TPS for the tissue by 
solving the matrix equation of Q = Ax for the vector x. 

11. A method in accordance with Claim 1 wherein the TBF is a 
cerebral blood flow (CBF). the TBV is a cerebral blood volume (CBV), the TMTT is 
a cerebral mean transit lime (CMTT), the TPS is a cerebral TPS, said determining a 
tissue type based on the TBF, the TBV, the TMTT, and the TPS comprises 
determining one of a viable tissue and a non-viable tissue based on the CBF, the CBV, 
the CMTT, and the cerebral TBS. 

12. A system comprising at least one of a computed tomography 
system and a nuclear magnetic resonance system, said system configured to: 

quantitatively determine a tissue blood flow (TBF) by deconvolving 
Q(t) and C a (t), where Q(i) represents a curve of specific mass of contrast, and C a (t) 
represents an arterial curve of contrast concentration; 
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quantitatively determine a tissue blood volume (TBV) by 
deconvoluting Q(t) and C a (t); 

quantitatively detemiine a tissue mean transit time (TMTT) by 
deconvoluting Q(t) and C £ (t); 

quantitatively detemiine a tissue capillary permeability surface area 
product (TPS) by deconvoluting Q(t) and C a (t); and 

determine a tissue type based on the TBF, the TBV, the TMTT, and the 

TPS. 

13. A system according to Claim 12 further configured to: 

quantitatively determine a TBV for a tissue having a blood stream 
containing a contrast without leaking the contrast into an interstitial space of the tissue 
by solving a matrix equation of Q = Ah for a vector h, wherein vector h includes a 
plurality of elements comprising an impulse residue function at different times, Q 
comprises a vector including elements comprising values of a tissue residue function 
at different times, A comprises a matrix formed by values of the arterial curve of 
contrast concentration at different times; 

quantitatively detemiine a TBF for the tissue by solving the matrix 
equation of Q = All for the vector h; and 

quantitatively detemiine a TMTT for the tissue by solving the matrix 
equation of Q = Ah for the vector h. 

14. A system according to Claim 13 further configured to 
determine a least squares solution for the vector h under an equality constraint. 

15. A system according to Claim 14 further configured to 
determine a least squares solution for the vector h under a time causality constraint 
and a minimum transit time constraint. 
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16. A system according to Claim 15 further configured to 
determine a hast squares solution for the vector h under a smoothness constraint, a 
nonnegativi: constraint, and a monotonicity ccmstraint, wherein the smoothness 
constraint foices h to be smoothly varying, and the monotonicity and nonnegativity 
constraints force h to start at a maximum and then monotonia) :!y decreases towards a 
zero baseline. 

1 7, A system according to Claim 12 further configured to: 

quantitatively determine a TBV for a tissue l aving a blood stream 
containing a contrast with leaking the contrast into an interstitial space of the tissue by 
solving a m:.:rix equation of Q = A>; from the linearization of the tissue residue 
function for ; vector x, wherein vector x includes a plurality of elements comprising 
TBF, TBV, "MTT, TPS and combinations thereof, Q comprises a vector including 
elements co::v:rising values of a tissue residue function ;it different times, A 
comprises a : itrix formed by values of the arterial curve o r contrast concentration 
and tissue r?: V.ue function at di :'.ererit times and combination? hereof; 

quantitatively dc: rminc a TBF for the tis: ; y solving the matrix 
equation of < 5 == Ax for the vector x: 

nantitatively determine a TMTT for the ti-nc 'n» solving the matrix 
equation of • • - Ax for the vector x; and 

i;unnti::!!ively d. rminc a TPS for the tiss y solving the matrix 
equation o r < - A:-: f-.r the ver x. 

IS. A system according to Claim 17 i.irther configured to 
determine a • t squares solut^n for the TBF, the TBV. - TMTT, and the TPS 
under a no-- '.ivitv constrain- that (determines a TBV, a V!?: ; : a TMTT and a TPS 
for the tissu' nere (here is le.rage of the contrast from ' lood stream into the 
interstitial s; 
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19. A system according to Claim 12 further configured to 
quantitatively determine a partial volume averaging scaling factor for the arterial 
curve of contrast concentration by: 

deconvolving the measured arterial curve of contrast concentration 
with a venous curve of contrast concentration to determine a transit time spectrum 
through a tissue of interest; 

extrapolating the arterial curve; and 

convolving the extrapolated arterial curve with the transit time 
spectrum to generate an extrapolated venous curve, wherein the partial volume 
averaging scaling factor is the ratio of an area underneath the extrapolated arterial 
curve to an area underneath the extrapolated venous curve. 

20. A system according to Claim 12, wherein the TBF is a cerebral 
blood flow (CBF), the TBV is a cerebral blood volume (CBV), the TMTT is a 
cerebral mean transit time (CMTT), the TPS is a cerebral TPS, said system further 
configured to determine one of a viable tissue and a non-viable tissue based on the 
CBF, the CBV, the CMTT, and the cerebral TBS. 

21. A computer readable medium encoded with a program 
executable by a computer for processing scanning data, said program configured to 
instruct the computer to: 

quantitatively determine a tissue blood flow (TBF) by deconvoluting 
Q(t) and C a (t), where Q(t) represents a curve of specific mass of contrast, and C a (t) 
represents an arterial curve of contrast concentration; 

quantitatively determine a tissue blood volume (TBV) by 
deconvoluting Q(t) and C a (t): 

quantitatively determine a tissue mean transit time (TMTT) by 
deconvoluting Q(t) and C a (t); 
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quantitatively determine a tissue capillary permeability surface area 
product (TPS) by deconvolving Q(t) and C a (t); and 

determine a tissue type based on the TBF, the TBV, the TMTT, and the 

TPS. 

22. A computer readable medium in accordance with Claim 21 
wherein said program is further configured to: 

quantitatively determine a TBV for a tissue having a blood stream 
containing a contrast without leaking the contrast into an interstitial space of the tissue 
by solving a matrix equation of Q = All for a vector h, wherein vector h includes a 
plurality of elements comprising an impulse residue function at different times, Q 
comprises a vector including elements comprising values of a tissue residue function 
at different times, A comprises a matrix formed by values of the arterial curve of 
contrast concentration at different times; 

quantitatively determine a TBF for the tissue by solving the matrix 
equation of Q = All for the vector h; and 

quantitatively determine a TMTT for the tissue by solving the matrix 
equation of Q = All for the vector h. 

23. A computer readable medium in accordance with Claim 22 
wherein said program is further configured to determine a least squares solution for 
the vector h under an equality constraint. 

24. A computer readable medium in accordance with Claim 23 
wherein said program is further configured to determine a least squares solution for 
the vector h under a time causality constraint and a minimum transit time constraint. 

25. A computer readable medium in accordance with Claim 24 
wherein said program is further configured to determine a least squares solution for 
the vector h under a smoothness constraint, a nonnegativity constraint, and a 
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monotonicity constraint, wherein the smoothness constraint forces h to be smoothly 
varying, and the monotonicity and nonnegativity constraints force h to start at a 
maximum and then monotonically decreases towards a zero baseline. 

26. A computer readable medium in accordance with Claim 21 
wherein said program is further configured to: 

quantitatively determine a TBV for a tissue having a blood stream 
containing a contrast with leaking the contrast into an interstitial space of the tissue by 
solving a matrix equation of Q = Ax from the linearization of the tissue residue 
function for a vector x, wherein vector x includes a plurality of elements comprising 
TBF, TBV, TMTT, TPS and combinations thereof, Q comprises a vector including 
elements comprising values of a tissue residue function at different times, A 
comprises a matrix formed by values of the arterial curve of contrast concentration 
and tissue residue function at different times and combinations thereof; 

quantitatively determine a TBF for the tissue by solving the matrix 
equation of Q = Ax for the vector x; 

quantitatively determine a TMTT for the tissue by solving the matrix 
equation of Q = Ax for the vector x; and 

quantitatively determine a TPS for the tissue by solving the matrix 
equation of Q = Ax for the vector x. 



27. A computer readable medium in accordance with Claim 26 
wherein said program is further configured to determine a least squares solution for 
the TBF, the TBV, the TMTT, and the TPS under a nonnegativity constraint that 
determines a TBV, a TBF, a TMTT and a TPS for the tissue where there is leakage of 
the contrast from the blood stream into the interstitial space. 
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28. A computer readable medium in accordance with Claim 21 
wherein to quantitatively determine a partial volume averaging scaling factor for the 
arterial curve of contrast concentration said program is further configured to: 

deconvolute the measured arterial curve of contrast concentration with 
a venous curve of contrast concentration to determine a transit time spectrum through 
a tissue of interest; 

extrapolate the arterial curve; and 

convolve the extrapolated arterial curve with the transit time spectrum 
to generate an extrapolated venous curve, wherein the partial volume averaging 
scaling factor is the ratio of an area underneath the extrapolated arterial curve to an 
area underneath the extrapolated venous curve. 

29. A computer readable medium in accordance with Claim 21 
wherein the TBF is a cerebral blood flow (CBF), the TBV is a cerebral blood volume 
(CBV), the TMTT is a cerebral mean transit time (CMTT), the TPS is a cerebral TPS, 
said program is further configured to determine one of a viable tissue and a non-viable 
tissue based on the CBF, the CBV, the CMTT, and the cerebral TBS. 

30. A computer readable medium in accordance with Claim 21 
wherein said program is further configured to: 

scan with at least one of a computed tomography (CT) system and 
Nuclear Magnetic Resonance system (NMR) to measure the arterial curve of contrast 
concentration and a tissue residue function; 

quantitatively determine a TBV for a tissue having a blood stream 
containing a contrast without leaking the contrast into an interstitial space of the tissue 
by solving a matrix equation of Q = Ah for a vector h, wherein vector h includes a 
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plurality of elements comprising an impulse residue function at different times, Q 
comprises a vector including elements comprising values of a tissue residue function 
at different times, A comprises a matrix formed by values of the arterial curve of 
contrast concentration at different times; 

quantitatively determine a TBF for the tissue by solving the matrix 
equation of Q = Ah for the vector h; 

quantitatively determine a TMTT for the tissue by solving the matrix 
equation of Q = Ah for the vector h; and 

quantitatively determine a TPS for the tissue by solving the matrix 
equation of Q = Ah for the vector h. 



31. A computer readable medium in accordance with Claim 21 
wherein said program is further configured to: 

scan with at least one of a computed tomography (CT) system and 
Nuclear Magnetic Resonance system (NMR) to measure the arterial curve of contrast 
concentration and a tissue residue function; 

quantitatively determine a TBV for a tissue having a blood stream 
containing a contrast with leaking the contrast into an interstitial space of the tissue by 
solving a matrix equation of Q = Ax from the linearization of the tissue residue 
function for a vector x, wherein vector x includes a plurality of elements comprising 
TBF, TBV, TMTT, TPS and combinations thereof, Q comprises a vector including 
elements comprising values of a tissue residue function at different times, A 
comprises a matrix formed by values of the arterial curve, of contrast concentration 
and tissue residue function at different times and combinations thereof; 

quantitatively determine a TBF for the tissue by solving the matrix 
equation of Q = Ax for the vector x; 
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quantitatively determine a TMTT for the tissue by solving the matrix 
equation of Q = Ax for the vector x; and 

quantitatively determine a TPS for the tissue by solving the matrix 
equation of Q = Ax for the vector x. 
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Rule 39.1(1v) PCT - Method for treatment of the human or animal body by 
surgery 

Rule 39 . 1 ( i v ) PCT - Diagnostic method practised on the human or 
animal body 
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because they are dependent claims and are not drafted in accordance with the second and third sentences of Rule 6.4(a). 
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